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Abstract 

We study the inclusive production of the next-to-hghtest neutralino X2' decaying 
directly, or via slepton, into two leptons and the lightest neutralino. The dilepton 
invariant mass spectrum in these decays has a characteristic sharp edge near the 
kinematical upper limit. We propose to exploit this feature as a search strategy 
for the X2' ^^^ thereby for SUSY. The possibilities to determine neutralinos and 
slepton masses are also discussed. 
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^ ' 1 Introduction 

Oh 



The Standard Model (SM), despite its phenomenological successes, is most 
likely a low energy effective theory of matter fermions interacting via gauge 
bosons. A good candidate for the new physics is the Supersymmetry (SUSY), 
which in its minimal version introduces a scalar (fermion) partner to each 
5^ \ ordinary fermion (boson) with the same couplings. If i?-parity is conserved, 

SUSY particles are produced in pairs and a stable Lightest Supersymmetric 
Particle (LSP) appears at the end of each sparticle decay chain. Weakly inter- 
acting LSP escapes detection and masses of sparticles cannot be reconstructed 
explicitly. Usually, one characterises the SUSY signal significance by an excess 
of events over the SM background expectation. 

SUSY partners of gauge and Higgs bosons mix to form physical states, charginos 
Xi,2i and neutralinos, Xi2 3 4- Depending on a particular SUSY scenario the 
decays X2~^ ^^^~Xi and/or X2~^ ^^^^ ~^ ^^^~Xi can take sizable branching ra- 
tios. The dilepton invariant mass spectrum in these decays has a specific shape 
with the sharp edge near the kinematical upper limit. In case of three-body 
(direct) decays X2~~* ^^^~Xi the position of edge is 

M12P = M^o - M^o (1) 
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and in case of two-body (cascade) decays it is 



^(M|o-Mf)(Mf-M|o) 

M,Tf- = ^ ^ (2) 

i+i ^ V ; 



As discussed in Ref. [1] for xf X2 Electro- Weak (EW) production, the three- 
body leptonic decays of X2 can be used to measure the mass difference be- 
tween X2 s-'^d Xi- At LHC, however, much higher production cross-sections 
are expected for strongly interacting gluinos and squarks. If q {g) is lighter 
than g (g), it can decay to quark plus neutralino or chargino: q—^ Qx'i/Q'xf 
{g^ qq'x'i/QQ'xf)- Eventually, this source of X2 is much more prolific than 
direct EW production. Therefore, we propose to perform inclusive search 
for X2 ^^^ ^'^ ^S6 the specific shape of the dilepton mass spectrum as evi- 
dence for SUSY [2]. For the SM background suppression, besides two same- 
flavour opposite-sign leptons, one can ask an additional signature character- 
ising SUSY. This can be: missing energy taken away by escaping LSPs, or 
additional jets coming from gluinos and squarks cascade decays, or an ex- 
tra high-j9T lepton from charginos, neutralinos, sleptons, IVBs, or b-quarks 
copiously produced in SUSY. Depending on sparticle masses and predomi- 
nant decay modes one of these extra requirements can turn out to be more 
advantageous than the others, or to be complementary. 

In this paper we discuss the possibility to observe the Xo, ^^ the inclusive 
3 leptons and 3 leptons + i?™*'' final states with the CMS detector at LHC. 
The search strategy for SUSY is to look for a characteristic shape in the 
dilepton invariant mass spectrum with a sharp edge, rather than just for an 
event excess over the SM expectations. Such an observation would reveal SUSY 
through X2 production, and at the same time would allow to constrain sparticle 
masses and some of the SUSY model parameters. We also propose a way to 
resolve the X2 decay type ambiguity (three-body versus two-body), and a 
method of slepton mass determination in the X2 two-body decay chain. 

The present analysis is performed within the framework of the minimal Su- 
persymmetry Model motivated by Supergravity (mSUGRA), where only five 
extra parameters are introduced [3]: m^, the universal scalar mass, nii/2, the 
universal gaugino mass and Aq, the universal trilinear term at GUT scale; 
tan(3, the ratio of vacuum expectation values of two Higgs fields; sign{fi), 
the sign of the Higgsino mass parameter. Once these parameters are specified 
all sparticle masses and couplings at EW scale are evolved via the renor- 
malization group equations. In the following we limit ourselves to the set: 
tanP = 2, Aq = 0, /i < and investigate what is the region of (mo, "021/2) 
parameter space, where the dilepton invariant mass edge is visible [2]. 



2 xi production and decay 

Within mSUGRA the following relations are valid: 

M^o ^ M^± ^ 2M^o ^ (0.25 ^ 0.35)M3 ^ 0.9mi/2 (3) 

M^ = m^ + 0.52 m2/2 - 0.5(1 - 2sm'^9w) m\ cos2/5 (4) 

Mf = m^ + 0.15 m^/2 " sin^^p^/ m| cos2/3 (5) 

The three-body leptonic decays of X2 dominate below rnxji ~ 200 GeV nearly 
for all values of mo, and there is also a small region at higher mi/2 where 
these decays are open. In these regions the measurement of MJ+f- within 
mSUGRA allows to get estimate of M^o (x?= LSP), M^o, M^±, Mg and 

mi/2 (c.f. eqs. (1), (3)). At mi/2 ~ 200 GeV, the modes, X2 -^ ^^Xi/^Xi, 
become kinematically allowed thus suppressing direct leptonic decays of X2- 
The two-body decays of the X2 ^^^ ^l/r occur mostly at a low values of rriQ. 
In these regions the measurement of only the edge position does not provide 
information on masses unambiguously (c.f. eq. (2)). However, as shown later, 
by analysing some kinematical distributions, the masses of Xi, X2 ^'^d of the 
intermadiate slepton can be determined with a reasonable precision. Therefore, 
with two-body decays one can constrain mi/2 and itlq if the information on 
tan/5 is available, e.g., from the Higgs sector of the model. In some regions of 
parameter space the two-body and three-body leptonic decays coexist or both 
modes X2 ~^ ^"^^r ^''^^ X2 ~^ ^^^l ^-re open, leading to a possibile observation 
of a " double edge" . 

At LHC the X2 i^ abundantly produced either from gluinos or from squarks 
over almost the whole {1710,1711/2) plane, as illustrated in Figs. la, 6. Figure 
Ic shows the X2 decay branching ratio into leptons (/ = /i or e). Finally, 
the X2 inclusive production cross-section times branching ratio into leptons 
is shown in Fig.lrf. Clearly, one expects quite large rates of opposite-sign 
same-flavour dileptons from inclusive X2 production over a large portion of 
(mo, mi/2)- Similarly to the X2, the lightest chargino is also copiously produced 
from strongly interacting sparticles. Moreover, the leptonic decay branching 
of the xt always exceeds 0.1 per lepton flavour. For low values of tuq, in the 
regions where decays to sleptons are kinematically allowed, the xf gives a 
lepton with a probability close to 1. Making use of this proliflc production 
of leptons from Xi , as well as from other SUSY sources, the SM background 
can be suppressed by asking for an additional high-p^ lepton. Thus the X2 
inclusive production can be studied in 3 lepton and 3 lepton + i?™*'' flnal 
states. 



mSUGRA parameters: tan[3 = 2, A^ = 0, |i< 
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Fig. 1. X2 inclusive production and decay: a) branching ratio 8(^^X2+^)) b) 
B{ul -^ X2+ '^)i ^) B(x2 -^ l^l^ + invisible) and d) X2 inclusive production 
cross-section times branching ratio into l^l^ (I = fi or e). 

3 Event simulation and selection 



For signal simulation we have generated all mSUGRA processes, using ISAJET 
7.14 [4]. The considered SM sources of background, WZ, ti, TIL and Z66 pro- 
ductions, have been simulated by PYTHIA 5.7 [5]. The following production 
cross-sections have been assumed at LHC energy of y^ = 14 TeV: (Tvfz = 26 
pb, Ott = 670 pb for a top mass of 170 GeV, azz = 15 pb and azbi = 580 pb. 



We have used a parameterised detector simulation program [6]. It takes into 



account: i) particle bending in a 4 T magnetic field; ii) smearing of lepton 
momentum according to parameterisation obtained from detailed simulations; 
iii) 90% triggering plus reconstruction efficiency per lepton within geometrical 
acceptance |?7| < 2.4; iv) 90% reconstruction efficiency per charged track with 
Pt > 1 GeV within |?7| < 2.4; v) coverage, granularity, main cracks, energy 
resolution and electronic noise of the calorimetry system; vi) fluctuation of 
starting point and the spatial development of electromagnetic and hadronic 
showers by parameterisation of the lateral and longitudinal profiles. For the 
high luminosity study, L = 10^'* cm~^s~^, on top of each signal and background 
events we have superimposed on average 15 "hard" pile-up events (PYTHIA 
QCD jet production with pT > 5 GeV) fluctuated by Poissonian distribution. 

In the inclusive 3 lepton channel we require: 1) two opposite-sign, same-flavour 
leptons (/i or e) with px > 10 GeV; 2) a third "tagging" lepton with px > 15 
GeV; 3) lepton isolation: if there is no track with px > 2 (1.5) GeV within 
AR = y^Arj'^ + A0^ < 0.3 around the lepton direction, it is considered as 
isolated. The leptons pr thresholds are variable depending on the region of 
(mo, "021/2) studied. In some cases leptons are not required to be isolated. In 
the inclusive 3 lepton + i?^*^*channel, beside the above three requirements 
we also ask: 4) EJp'''> 200 (300) GeV. 

In the selected events we reconstruct the invariant mass Mi+i- of the lepton 
pair with the same flavour and opposite sign. When several /"''/" combinations 
per event are possible (in case of, e.g., three leptons of the same flavour) the 
one with the minimal distance Ai?^ ' is chosen. 



4 Dilepton invariant mass spectrum 



The dilepton invariant mass spectrum for a representative mSUGRA param- 
eter space point {niQ = 200 GeV, nii/2 = 100 GeV) superimposed on the ex- 
pected SM background is shown in Fig. 2a. The number of events corresponds 
to an integrated luminosity of Lint = 20 pb^^, i.e. a few hours of LHC running 
at a luminosity of L = 10^'^ cm^^s^^. The pt thresholds on all three leptons 
are 15 GeV and two of them entering mass distribution are isolated. The SM 
background is negligible, with the dominant contribution from WZ produc- 
tion around the Z mass peak and ti production at lower invariant masses. The 
specific shape of the distribution with its sharp edge reveals SUSY through X2 
production. At this mSUGRA point the X2 ^as three-body decay modes, and 
the relevant masses are M^o = 97 GeV, M^o = 45 GeV. The edge is situated 
at the expected value of MJ+I- = M^o — M^o = 52 GeV and its position can 
be measured with an accuracy better than 0.5 GeV. Within the model this 
allows to determine M^o, M^o, M-±, Mg and mi/2 (c.f. (3)). Note, that the 
Z peak seen in Fig. 2 serves as an overall calibration signal; it allows to control 
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Fig. 2. Expected I'^l^ mass spectrum for mSUGRA points: a) niQ = 200 GeV, 
777-1/2 = 100 GeV in the inclusive 3 lepton final states; b) mo = 150 GeV, 7771/2 = 700 
GeV in the 3 lepton + £^™** final states. The other mSUGRA parameters are 
tan/3=2, ^0 = and ^ < 0. The hatched histogram corresponds to the SM back- 
ground. 

the mass scale as well as the production cross-section. 

With increasing 7721/2 the observation of the "edge" becomes more difficult 
due to a rapid fall of gluino/squark production cross-sections. Here, on the 
other hand, heavier g and q provide higher missing transverse energy. Asking 
EJp^^^> 200 (300) GeV rejects most of the SM background whilst retaining 
a large fraction of signal events. Figure 26 shows the dilepton spectrum for 
the mSUGRA point (150 GeV, 700 GeV) in the 3 lepton + E^''' final states. 
In this region of high mi/2 dileptons are mainly produced through cascade 
decays of X2 through intermediate state sleptons. The third, "tagging" lep- 
ton, predominantly comes from a Xi produced in association, again from g, q 
cascades. 



In order to delineate the region of parameter space where the /"'"/" edge is vis- 
ible, the {1710,1711/2) plane has been systematically scanned for fixed tan/? = 2, 
Ao = and /i < 0. At least 40 signal events in a mass interval > 10 GeV just 
below the edge position and a statistical significance S = Ns/^/Ns + Nb > 7 
are required {Ns and A^^ are the number of signal and background events in 
this mass interval). Figure 3 shows the domain explorable under these condi- 
tions with an integrated luminosities of Lint = 10^ pb~^ and Lint = 10^ pb^^ 
for 3 lepton and 3 lepton + EJp'^^'^ channels, respectively. Of particular impor- 
tance is the fact, that already at Lint = 10*^ pb~^ the cosmologically preferred 
region 0.15 < Vth? < 0.4 [7], where the x? would be a good candidate for dark 
matter, is entirely covered; here Vt is the ratio of the x5 "^^^^c density to the 



mSUGRA parameters: tan^=2, Ag=0, \x<0 
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Fig. 3. Explorable region of niSUGRA in inclusive 3 lepton and 3 lepton + E"™**^ 
final states with Lint = 10^ pb^^ and Lint = 10^ pb~^, respectively. Shaded regions 
are excluded by theory and/or experiment. The simulated mSUGRA points are 
shown by circles. The cosmologically preferred region 0.15 < il/i^ < 0.4 is also 
given. 

critical density of the Universe and h is the Hubble constant scaling factor. In 
the 3 lepton + EJp^'^^ final states with Lint = 10^ pb~^ the reach extends up to 
mi/2 ~ 900. A detectable edge is seen as long as a-B > 10~^ pb (see Fig. Id). 



5 X2 decay type and slepton mass determination 



The X2 decay type (three-body versus two-body) determination is essential 
for the interpretation of the observed edge. The px spectra of leptons from 
X2 (cascade) two-body decays are more asymmetric compared to the three- 
body ones. To characterise this asymmetry we introduce the variable A = 
where p™"^ (p™*") corresponds to the lepton of maximal (minimal) 



Pt 
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transverse momentum. An example of decay-type determination is given for 
the mSUGRA point (100 GeV, 150 GeV). Here the X2 decays via X2 ^ l^^R ^ 
l~^l~X^ with 0.54 branching ratio and the sparticle masses are M^o = 135 GeV, 



M.o = 65 GeV, 



Mr 



120 GeV. Figure 4 shows the dilepton spectrum for 
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Fig. 4. I'^l mass distribution for mo Fig. 5. Lepton pj'-asymnietry distribu- 



= 100 GeV, mi/2 = 150 GeV, tan/3=2, 
^0 = s-nd // < 0. The hatched his- 
togram corresponds to the SM back- 
ground. 



tions in three- and two-body decays of 
X2 with the same position of dilepton 
invariant mass edge. 



this point with Lint = 5 ■ 10^ pb^^, superimposed on the SM background, 
in the inclusive 3 lepton final states. An edge is situated at ~ 52 GeV, as 
expected. To look at the lepton p-r- asymmetry distribution, we pick up the 
lepton pairs with Mi+i- < 52 GeV. Figure 5 (full line) shows the corresponding 
A spectrum for the selected events. The pronounced asymmetry in px indicates 
the cascade nature of X2 decays. The dotted histograms are obtained with the 
assumption that the "observed" edge is due to direct three-body decays with 
various M^o, M^o combinations providing M^o — M^o = 52 GeV. In these cases 
the asymmetry A distributions peak at zero. 



After the decay type is determined a further analysis step is carried out to 
extract the masses of involved particles: i) assume M^o = 2M^o, ii) generate 
samples of X2 2-body cascade decays for various M^o (M^); note, that the 
slepton mass is constrained to provide the "observed" position of an edge 
given by eq. (2), iii) the best combination of M^o, M^o and Mj is then chosen 
by means of a x^-test of the shape of the lepton pj^- asymmetry distributions. 
This procedure provides the following precisions on masses: SM^o < 5 GeV, 
6Mj < 10 GeV. The use of the total number of observed events, as well as the 
dilepton invariant mass spectrum itself in a combined x^-test would further 
improve these results. 
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Fig. 6. Predicted l^l^ mass spectra for 
two values of slepton mass: Mj = 120 
GeV (full line) and Mj = 77 GeV 
(dashed line) ; mSUGRA parameters are 
as in Fig. 4. 
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6 Double edge 



As discussed in section 2, in some regions of the (mo, fni/2) plane the different 
leptonic decay modes of X2 can be simultaneously open and observation of 
"double edges" is possible if the branching ratios are not too dissimilar. This is 
the case, for instance, for the mSUGRA point (100 GeV, 150 GeV) discussed in 
the previous section. Here, beside the X2 ~^ ^^^r ~^ ^^^ X? decays, the direct 
three-body decay mode is also open, though, with a much smaller branching 
ratio, B{xl -^ l^l'Xi) = 0.05. The second edge at M^'^ff = 69 GeV seen 
in Fig.4 is thus less spectacular than the one at 52 GeV. The px-asymmetry 
distribution of dileptons with 55 GeV < Mi+i- < 69 GeV indicates the three- 
body decay type of these events. Using now two measured values of edge 
positions and assuming M^o = 2M^o, we obtain two solutions for the slepton 
mass from equations (1) and (2): Mj = 120 GeV and 77 GeV. Corresponding 
dilepton invariant mass spectra for these two solutions are shown in Fig. 6. The 
clear difference in the predicted spectra (Fig.6 vs. Fig.4) allows to eliminate 
the Mj = 77 GeV solution. At this particular mSUGRA point the use of both 
edges provides a precision of 6M-0 [ ^ 1 GeV. 



Another example of a "double edge" is given in Fig. 7 for the mSUGRA point 
(50 GeV, 125 GeV). Here the sparticle masses are M^p = 116 GeV, M, 



'xi 



55 



GeV, Ml = 110 GeV, Mj = 78 GeV and now two two-body decays coexist, 
with a comparable branching ratios B(x2 -^ ^^^l ~^ ^~^^~Xi) = 0.037 and 
B(x2 -^ f- Ir -^ ^■~^J-~Xi) — 0.013 and their analysis could proceed as indicated 



above providing a strong constraint on the underlying model. 



7 Conclusions 



• Observation of an edge in the dilepton invariant mass spectrum reflects 
production of X2 ^^^ hence establishes existence of SUSY; this observation in 
some cases is possible with very small statistics and could be the flrst evidence 
for SUSY at LHC. 

• In the inclusive 3 lepton flnal states, with no great demand on detector 
performance, a large portion of parameter space, including the cosmologically 
preferred mSUGRA domain, can be investigated at an integrated luminosity 
of Lint = 10^ pb^i. 

• In the inclusive 3 lepton + E"^*^* flnal states and with Li^t = 10^ pb~^ most 
of the region where the X2 inclusive production cross-section times branching 
ratio into leptons > 10~^ pb can be covered. 

• From the observation of characteristic edge(s) and with mSUGRA motivated 
assumptions, it is possible to determine masses of Xi and X2 i^i case of direct 
decays and masses of Xi, X2 ^^^ sleptons in case of cascade decays. 
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